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Abstract 

The free fermion condition of the six-vertex model provides a 5 parameter sub- 
manifold on which the Bethe Ansatz equations for the wavenumbers that enter into 
the eigenfunctions of the transfer matrices of the model decouple, hence allowing 
explicit solutions. Such conditions arose originally in early field-theoretic S-matrix 
approaches. Here we provide a combinatorial explanation for the condition in terms 
of a generalised Gessel-Viennot involution. By doing so we extend the use of the 
Gessel-Viennot theorem, originally devised for non-intersecting walks only, to a 
special weighted type of intersecting walk, and hence express the partition function 
of N such walks starting and finishing at fixed endpoints in terms of the single walk 
partition functions. 
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1 Introduction 



There has been a steady stream of interest in the statistical mechanics of directed walk 
problems because of their connections to the physics of polymers and domain walls |T|. 
Since the popularisation of the field in the seminal article by Fisher vicious walkers, 
in particular, also known as non-intersecting walks, on two-dimensional (directed) lattices 
have been the subject of much work p, ^ |[ |], |], 0]. In the field of combinatorics a very 
general methodology, useful for any directed graph and based on a involution, has been 



developed by Gessel and Viennot P, H following the work of Lindstrom []TU| , and Karlin 



and McGregor [|lT|, which expresses the generating function of configurations of walks 
as the value of a determinant of single walk generating functions. In the most general 
setup an arbitrary inhomogeneous weight may be associated with each occupied edge of 
the lattice. 

On the other hand it has been well known for a long time ||12|, |^ that the square 
lattice six-vertex model can be mapped onto a problem of interacting (intersecting) di- 
rected walks on that lattice (see figure m). Let us call the six weights of that model 
{wi,W2,W2,,Wi,W5,WQ}: see figure |1]. We can consider wi = 1 without loss of generality. 
In order to calculate the partition function of fixed numbers of walks one needs to consider 
a particular invariant sub-space of the associated transfer matrix, the diagonalisation of 
which involves the famous Bethe Ansatz trial solution: The Bethe Ansatz is a guess for 
the eigenvectors of the transfer matrix and involves a sum over a set of plane wave forms. 
For a problem of walks the Bethe Ansatz involves wavenumbers which are chosen 
from the solutions of a set of non-linear coupled polynomial equations. To calculate 
the walk partition function one needs to find, and be able to sum over, all the eigenvalues 
and eigenvectors explicitly. 

Recently it has been shown that if one rather considers the combinatorialist's 



problem of A^ vicious walkers with weights associated with edges, rather than with ver- 
tices, this can be solved using the transfer matrix/Bethe Ansatz approach in a completely 
rigorous fashion. Here the Bethe Ansatz equations for A^ walks decouple and the solution 
of the A^ walk problem is given by the Gessel- Viennot determinant of single walk gener- 
ating functions. The walk problem using edge weights is equivalent to a restricted vertex 
model (or visa versa) where the vertex weights must satisfy the equations 

WsWi = W^Wq (1.1) 

and 

W2 = (1.2) 
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The second condition ( |1.2|) merely expresses the fact that the weight associated with the 
meeting of two walks is set to zero since vicious walkers are being considered. The asso- 
ciated vertex model with only this second condition necessarily holding is often referred 
to as the five vertex model. 

Of central importance here is that taken together the two conditions, (|1 . 1| ) and ( |1.2|) , 
imply that the less restrictive free-fermion condition of the six-vertex model, which occurs 
when the vertex weights satisfy the equation 

W1W2 = - W3W4, (1.3) 

is then automatically satisfied. Hence the edge weight model is equivalent to the free- 
fermion case of the five-vertex model. This is not surprising since the free-fermion condi- 
tion is precisely the general condition needed to achieve the decoupling of the Bethe 
Ansatz equations in the solution of the six-vertex model. 

This raises the question of whether the Gessel-Viennot methodology can be adapted to 
the 'free-fermion' case (i.e. equation ( |1.3| ) is satisfied), of 'six- vertex' or 'osculating' walks 
- here the walks are allowed to intersect but not share edges i.e. site-only intersecting. 




W-^ W2 Wq 

Figure 1: At the top are the six possible path configurations at a vertex of the lattice. 
Below each of these is the associated arrow configuration of the six-vertex model, while 
below that is the six weights we associate with each of those configurations. 

In this paper we demonstrate that free-fermion osculating walks can indeed be counted 
with a generalisation of the Gessel-Viennot methodology and we hence explicitly calculate 
their generating function. Because the Gessel-Viennot involution is involved the result is 
again a determinant of single walk generating functions. We restrict our discussion to the 
square lattice but the ideas can be easily generalised to any directed (acyclic) graph. In 
order to understand the combinatorial interpretation of the osculating free-fermion walks 
we briefiy discuss the case of non-intersecting walks for those unfamiliar with the Gessel- 
Viennot method. For the case of non-intersecting walks the Gessel-Viennot theorem gives 
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Figure 2: An example of three osculating paths. 

the generating function of such paths starting and ending at some fixed sets of sites as 
the determinant of the generating functions of independent-single-walks where the walks 
start and end at all permutations of those endpoints. That is, the generating function, 
■^t^\'y^ ^ 'y'^)^ for ^ non-intersecting walk configurations where the j*^ walk starts at 
i/j and arrives at after t steps is given by the following determinant: 

A/'f^(y^-yO = det||M«(y^^y;^)|| a,l3=l...N- (1-4) 

where = {y\, . . . , y]^) and y-^ = {y(, . . . , y^). Hence the generating function of walks 
is given as the signed sum over products of single walk generating functions where the 
sign of the contribution is signature of the permutation of the endpoints. This signed 
sum is interpreted as a signed sum over the elements of a "signed" set, Q, the elements 
of which are configurations of A^, possibly intersecting, walks, 

N 

^r(y^-yO = E^xn-^J'^(^")- (1-5) 

xen a=l 

where Xa is a single walk from Ua^yp, ex is the sign of the configuration and A4\ (Xa) 
the weight of the a particular configuration of a single walk - see section § for more 
precise definitions of the terms referred to in this section. Note that the A^ single walks in 
elements of Q may be edge as well as site intersecting, as they are completely independent 
of each other. 



The method introduced by Gessel and Viennot 0, M shows that pairs of opposite 



signed terms of (L5) corresponding to configurations with at least one intersection cancel 
leaving only positively signed terms corresponding to non-intersecting configurations - 
precisely the sum required to give J^t^\'y^ -^y-^)- This pairing rule is an "involution". 

In order to interpret the free-fermion condition for osculating or six-vertex walks we 
apply the same type of methodology but now require an extended pairing rule, one for 
which the pairs of terms in the signed sum do not necessarily cancel out. We now consider 
three different disjoint subsets of Q: Q"' the subset containing configurations that do not 
intersect at all, the subset containing configurations that only share sites (and not 
edges) and Q"^ the subset containing configurations that share at least one edge. 

The involution now pairs oppositely signed terms from Q'^ which cancel out, but the 
pairs of oppositely signed terms from the subset have different weights and do not 
cancel out, rather the weight difference is precisely given by the free-fermion condition 
( |1.3|) . This idea is expressed schematically by the case for (example) paths in the subset 



il'' as 

nX)-nX) = W(X) (1-6) 
and for walks in the subset fi^ as 



V{^)-V{y^) = (1.7) 

where V{X) is the product of the vertex weights of the two single paths of the configuration 
X, whilst VV(X) is the vertex weight of the paths, X, taken as a whole. 

Thus the free fermion condition arises as a natural consequence of summing over a 
signed set of A^-walk configurations each of whose weight is a product of single path vertex 
weights. The intersecting configurations in the signed set which do not cancel combine to 
create the correct vertex weight for the iV-walk configuration. 



2 Free Fermion walk generating function 

Definitions and Notations We will consider walks on the directed square lattice, C 
rotated through 45°. Each vertex of the lattice is labelled by the "time" coordinate, t 
and a height y and represented by the pair (t, y) or function y{t) - see figure |[ An 
N -vertex is an A^-tuple of vertices of C each of has the same time coordinate and is 
represented by the A^-tuple of their height coordinates, y = (yi, . . . , y^)- An A^- vertex y 
is non-intersecting ii yp ya, a, f3 = 1, . . . , N, a p. A walk of length t, Xiija^yp) = 
y{0)y{l)y{2) . . .y(t), where y{0) = y^ and y{t) = yp, is a sequence of t adjacent edges 
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from vertex y{0) to vertex y{t). An A^-walk, X(y*, y-^) is an A^-tuple of equal length walks, 
{Xi{y\,y(), . . . ,XN{yN,yN)), with and y-^ non-intersecting. 

An osculation between walks Xa and Xj3 at time s occurs if ya{s) = Visis), and 
either y^is - 1) < yfsis - 1) and ya{s + 1) < V/sis + 1) or y^is - 1) > yfsis - 1) and 
ya{s + 1) > yi3{s + 1) occur i.e the configuration xx appears. Paths X^ and X^ switch 
at time s if = yjsis), and either ya{s — 1) < yp{s — 1) and y^is + 1) > + 1) or 
ya{s — 1) > yfs{s — 1) and ya{s + 1) < + 1) occur i.e the configuration X appears. 
Since the initial and final A^-vertices are non-intersecting we do not define any osculations 
or switches associated with these vertices. A pair of walks is osculating if all common 
vertices are osculations - see figure 0. An A-walk X is osculating if all the vertices in 
common with any of the walks form osculations. 

Let Xa = ?/„(0)y„(l)?/a(2) . . . ya(t) and X^ = yp{0)yp{l)yp{2) . . . yp{t) intersect at time 
s i.e ya{s) = yp^s), then the pair X^ and Xp are exchanged at time s if 



Define [A^] = {!,... , A^}. Let P/v be the set of permutations of [A^], then for a = 
(ai, . . . , (TAr) e Pn, and y = {yi, ... , t/at) an A^-vertex, then cr(y) = {y„^,. . . , y„^). The 
signature of a permutation is denoted e^-. 

We will associate weights with the walks on the lattice as follows. Associate a set of 
six vertex weights, 



with each vertex v G C. Without loss of generality we only consider the situation where 
five of the weights are not equal to one. The vertex weight, W(X) of a particular A-walk, 
X, is the product of the vertex weights of all the vertices of C that are traversed by the 
walks of A^-walk. For each vertex traversed by the A-walk, only one of the six possible 
vertex weights associated with the vertex of C is used. Which of the six possible weight 
used depends on which of the four edges adjacent to the particular vertex are traversed 
by the A^-walk as illustrated in figure |l|. No vertex weights are associated with the initial 
and final A-vertices of the A^-walk. 

Remark. Note that with the above definition of osculating walks the set includes non- 
intersecting walks also. 

Remark. The vertex weight VV(X) is only defined for osculating paths., i.e if the paths 
of the A^-walk do not cross or if X contains only one path. 



Xa^ X'^ = ya{0) . . . ya{s)y(3{s + 1) . . . yfs{t) 
Xp-^ X'p = y/siO) . . . yf}{s)ya{s + 1) . . . ya{t) 



(2.8) 



V{v) = {wi{v), . . . ,we{v)}, 



(2.9) 



Remark. This way of associating weights with the walks is a generahsation of the sit- 
uation dealt with by the Gessel-Viennot theorem where the weights are associated with 
the edges of the graph and are independent of the A^-walk configuration. The problem 
of finding the generating function for osculating walks with arbitrary weights requires a 



rather complicated involution |Tj], however in the special case where the vertex weights 
satisfy the free fermion equation the osculating walk generating function can be expressed 
ClS db determinant, as we demonstrate here. 

Theorem [Free Fermion Walks]. Let Q* be the set of all osculating N -walks of length 
t starting at y' and ending at y-^ with yl^ < yl^_^_i and yj^ < yl^^i, a G [A^ — 1]. //, for 
V E C, the vertex weights satisfy 

W5{v)we{v) - W3{v)w4{v) = W2{v) (2.10) 

then the osculating lattice walk generating function, 

^fHy^-yO= E (2.11) 

where W(X) is the vertex weight of the N-walk X, is given by 

J^^ir - yO = Det\\^i'\yl - yi)\U=i...N (2.12) 

where J^t^\ya — > y^) is the generating function for a single lattice walk from y^^ y^- 

The theorem is proved by an extension of the Gessel-Viennot involution which does 
not preserve the A^-walk weight. 

Proof. Consider the signed set Vt = U , Q'^ fl = (p, where (p is the empty set. 
The positive and negative sets are 

n+ = {X(y'',a(yO)|a G Pn and = +1} (2.13) 
n- = {X(y^a(y■^))|^T G Pn and = -1} (2.14) 

Let, X G fi, then the sign of X is defined as 

f +1 if X G , , 

[ -1 if X G fi" 

For the A^-walks of Q we do not use the vertex weight of the A^-walk as a whole, but 
rather define a "product" weight, V(X.) in terms of its individual walks. In particular, 

P(X) = II W(X„) (2.16) 



We construct a sign reversing involution^, -0 on Q. 

Let X = {Xi, . . . , Xjv) G Q. The involution is an extension of the Gessel-Viennot 
involution and splits into three cases, 

1. No intersections (X e If none of the walks of X intersect, then ip(X.) = X. 

2. Edge intersections (X G Q'^). If any of the walks G X has an edge of C in 
common with another walk, X^ G X then let a be the least integer for which a walk 
Xa shares an edge with another walk Xf3. Of all the vertices adjacent to the edges in 
common with Xa and Xf3 choose the one with the smallest t coordinate and denote 
it by V, then X' = '0(X) is defined as the A^-walk obtained by exchanging X^ and 
Xp at V. 

3. Vertex only intersections (X G Q^). If any of the walks of X share vertices and none 
share edges, then let a be the least integer for which a walk X^ intersects another 
walk Xi3. Of all the vertices in common with Xa and X/3 choose the one with the 
smallest t coordinate and denote it by t;, then X' = ^/'(X) is defined as the A^-walk 
obtained by exchanging Xa and Xp at v. 

Call the vertex, v at which the involution exchange takes place, the "involution" vertex. 
The difference in the product weight of X and is then 

(w5{v)wq{v) — W3{v)w4{v)) P{^}(X) if V is an osculation 
{w3{v)w4,{v) — W5{v)wq{v)) V{.u}PQ if is a switch 
otherwise 



P(X)-W(X)) 



(2.17) 

where V^^y (X) is the product weight of X with the contribution of the weight associated 
with V divided out. 

Since C is defined as the set of A^-walks for which the involution vertex v exists 

and arises from the "vertex only intersections" case of ip, this means that all the A^-walks 

in have walks which only intersect at vertices - there are no shared edges. Define two 

A^-walks to be related, X ~ X' iff X can be obtained from X' by the interchange of any 

number of osculations with switches (or visa versa). This relation is easily seen to be an 

equivalence relation and hence partitions Q'^ into disjoint subsets, fi^, a G X where X is 

sign reversing involution, ■0 is a permutation of fl such that ip'^ = Identity and it has the property 
that whenever ■0(X) ^ X, then X G if and only if ipiX) e il^ . 
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some index set for the partitions. Define the canonical element, X'^'* of each partition, 
f2^, as the A^-walk for which all the intersections are osculations. Let S(X.'^") be the set 
of vertices of X*^" G fi^ in common with at least two walks of X^". Note, the cardinality 
of is2l^(x^")l. 

For the A^- walks in we have, 

5^ exP(X)=P^(X^'^) H {w,iz)we{z) - ws{z)w,{z)) (2.18) 

where 'P£-(X^") is the product weight of X'^" with all the vertex weights associated with 
the vertices in £{'K'^°') divided out. This follows since Y[z€£{X'=a){'^5{z)w6{z) —W3{z)w4{z)) 
allows for each vertex in S(X.'^°') to be a switch (i.e. weight W3W4) or an osculation (i.e. 
weight w^Wq). The sign, ex is correctly obtained since it is just —1 to the number of 
occurrences of a switch i.e. the number of factors of —w^w^. Thus we have the following: 



Det\\V{yl, ^ yi)\Up=i...N = Yl exP(X) (2.19) 

= Yl ^xP(X) + Y ^x^(X) (2.20) 

xen-n= xgf^s 

and since, by (|2.17|) the A^-walks in — (= U Q"^) with any intersections cancel in 
pairs we get 

= J2 ^(^) + E ^x^(X) (2.21) 

xefi" xgH'* 

= J2 + E E ^x^(X) (2.22) 

where C f2 is the set of non-intersecting A^-walks, using (|2.18|) gives, 

= J] >V(X) + ^P^(X'=") H {w,{z)w,{z) - ws{z)w,{z)) (2.23) 
xeH" aei 26£:(x=«) 

now, using the free fermion relation, (|2.10|) we get 

= 5^ W(X) + 5^P,(X-) n ^2(^) (2.24) 
xeH" aei zg£:(X':<^) 

= ^ >V(X) + ^W(X^") (2.25) 

xef^" aex 

and since all the canonical A^-walks, X^° are osculating we get 

= Y (2-26) 

xefi* 

as required. 

□ 
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